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QCD COLOR GLASS CONDENSATE MODEL IN WARPED BRANE MODELS 
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Hadron-hadron interaction and Deep Inelastic Scattering (DIS) at very high energies are dominated by events at 
small- xb regime. Interesting and complex physical content of this regime is described by a phenomenological model 
called McLerran-Venugopalan Color Glass Condensate (MVCGC) model. The advantage of this formalism is the 
existence of a renormalization-type equation which relates directly observable low energy (small- ) physics to high 
energy scales where one expects the appearance of phenomena beyond Standard model. After a brief argument about 
complexity of observations and their interpretation, we extend CGC to warped space-times with brane boundaries 
and show that in a hadron-hadron collision or DIS all the events - and not just hard processes - have an extended 
particle distribution in the bulk. In other word, particles living on the visible brane escape to the bulk. For an 
observer on the brane the phenomenon should appear as time decoherence in the outgoing particles or missing energy, 
depending on the time particles propagate in the bulk before coming back to the brane. Assuming that primaries of 
UHECRs are nucleons, the interaction of Ultra High Energy Cosmic Rays (UHECRs) in the terrestrial atmosphere 
is the most energetic hadron-hadron interaction available for observation. Using the prediction of CGC for gluon 
distribution as well as classical propagation of relativistic particles in the bulk, we constrain the parameter space of 
warped brane models. 
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1. Introduction 

One of the wondering fact about the Nature as we know 
it, is the huge difference between the strength of the 
Gravity and other forces. Another manifestation of this 
puzzle is called mass hierarchy - it can induce large ra- 
diative corrections to the mass of Standard Model (SM) 
.particles in its simplest Grand Unification (GUT) ex- 
tension. In recent years the idea of a TeV scale gravity 
in the context of large extra dimensions and localized 
'matter on the brane boundaries raised by N. Arkani- 
JHamed, et al. ^ and by L. Randall and R. Sundrum |2| 
inspired by some previous works by V. Rubakov and 
M. Shaposhinkov [3] on domain walls in higher dimen- 
sional spaces and by L Antoniadis |2,P. Horava and E. 
Witten |S] on M-theory models with compactification 
in spaces with D-brane boundaries, have created lots 
of excitements and hope for solving this long standing 
puzzle. 

In the first proposals it was suggested that only grav- 
ity should propagate in the higher dimensional bulk. 
Later it was proved that a total localization of all the 
fields except graviton on 3-branes is not realistic. In 
fact brane solutions are cosmologically unstable and at 
least one scalar bulk field (radion) [HI El is necessary 
to stabilize the distance between branes. In some brane 
models inflaton |S] also has to propagate to the bulk 
to make inflation with necessary properties. A deeper 
insight to the propagation of gravitational waves and 
massive particles with bulk modes in models with infi- 
nite bulk has illustrated that even the warping of the 



^e-mail: hzOmssl.ucI.ac.uk 



bulk can not stop their escape from branes |H] JOI- In 
fact geometric confinement in warped models is not as 
efficient as it was expected and specially spin-1 gauge 
fields jllj can propagate (tunnel) to a warped infinite 
or macroscopic bulk |EI ■ Due to violation of the gauge 
symmetry ,13) ,14), it is practically impossible to local- 
ize gauge field completely and only some modification in 
their interactions such as enhancement of their coupling 
to localized fermions on the visible brane can partially 
confine them ^2] QZI ■ Even this remedy is not with- 
out drawbacks and the universality of fermions' charge 
can be violated . Thus, it does not seem to be possi- 
ble to confine all the SM fields on the brane even when 
the bulk metric is warped geometrically or just by a mild 
modification of the Lagrangian. The only way out - if 
branes really exist - is a symmetry which keeps at least 
most fields at energies lower than a threshold confined 
to the branes. This would be possible if branes are topo- 
logical defects related to quantum gravity and created 
during a phase transition epoch in the early universe. If 
this broken symmetry prevents the production of KK- 
modes at low energies e.g. energies lower than the scale 
of broken symmetry, their direct production at present 
accelerator energies would be completely forbidden. For 
the same reason their effect on Zq/W^ width as well as 
gluon propagator would be extremely suppressed. For 
not having to add another scale to the theory, this scale 
should be close to A/5 , the fundamental scale of gravity. 

Propagation of particles in the bulk has a number of 
cosmological consequences which apriori can be used to 
constrain brane models. Nevertheless, the uncertainty 
of cosmological measurements and the dependence of 
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interpretations on the cosmological models does not yet 
permit to rule out these models completely. Thus a 
more controllable and/or close to home test is highly 
appreciated. 

Constraints on the parameter space of brane world 
models from collider data is mainly based on the prob- 
ability of direct observation of processes involving the 
production of gravitons and its Kaluza-Klein modes ^Hl 
[T^ [201 [21] ■ The accessible energy scales to the exis- 
tent and near future accelerators is however limited and 
the fundamental scale of gravity and the effective size 
of the extra-dimensions can be constraint only up to 
~ SOTeF. 

Another possibility is trying to find a signature of 
TeV scale gravity and large higher dimensions in the 
high energy air showers produced in the terrestrial atmo- 
sphere by Ultra high Energy Cosmic Rays (UHECRs). 
Assuming that the primaries are single elementary parti- 
cles (most probably proton or anti- proton) , these events 
are the most energetic particle collisions available to us 
at present. The energy scale of these events is ~ 10^^ eV 
roughly 3 orders of magnitude higher than available en- 
ergies even in near future accelerators such as LHC. 
If the fundamental scale of gravity in brane models is 
around Electroweak scale, UHECR showers with Cen- 
ter of Mass (CM) energies larger than > ITeV are 
expected to overcome symmetry restrictions and extra- 
dimension(s) should highly influence the behavior of 
showers. 

To better understand the potential of UHRCRs as a 
laboratory for observing new physics, we should remind 
ourselves that there is an essential difference between 
observables in colliders and in the air-showers. In the 
former case only particles with a transverse momentum 
greater than a minimum value which depends on the 
detector hole size are detectable and the remnants of 
the colliding beams which include most energetic par- 
ticles are not visible. In contrast, in an air-shower it 
is apriori possible to detect all the particles specially 
the most energetic ones and there is no discrimination 
between semi-hard and high energy remnants. Conse- 
quently, one is not restricted to see only high transverse 
momentum part of the collision. The remnant of the 
hadron after the collision consists of very energetic par- 
ticles which come from the scale to which the hadron was 
boosted. They can carry important information about 
these scales which are usually unobservable in labora- 
tory colliders. 

Treating particles classically, their propagation in 
bulk leads to a delay in their arrival if they are reflected 
back by the gravitational warping. When this time delay 
is very short, this phenomenon is equivalent to a larger 
effective mass (heavy Kaluza-Klein modes). If however 
the time delay is long, it appears as a time decoherence 
in the shower - some particles arrive detectably later 
than others. Due to longer fly distance of particles of 
the air showers produced by UHECRs, they are more 



sensitive to decoherence than colliders. 

If during propagation into the bulk particles are con- 
sidered as free, their wave equation can be solved and 
the spectrum of KK-modes can be determined. In the 
warped brane models the spectrum of KK-modes is very 
close to a continuum beginning from the 5-dim mass 
with a slight gap between the zero-mode and higher KK 
modes [H [12 [HI [ISl [23- The mass of KK-modes 
for all types of fields has the following general form 

[H [22 [13 [23 : 
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where m„ is the mass of n*'* KK-mode, a;„ > 3 for n > 
and its exact value depends on the spin of the field. 
Parameter L is the effective size of the bulk. The gap 
between the zero-mode and higher modes is determined 
by the scale of the compactification /j, . Warping of the 
extra-dimension is essential for explaining the observed 
hierarchy between electroweak scale and 4-dim Planck 
scale. In the RS-metric, warping is the conformal factor 
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ds^ = —irjn.vdx^dx'^ — dz^) 
e-^^'y^]f,^dx^'dx'' - dy^ 
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where branes are considered to be at z = R — l/fj, and 
z = R' = —e'^^ , and L is the distance between two 
branes. Warping leads to a mass spectrum for RS-type 
models where the KK-modes can be quite light: 
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If we consider /i as a free parameter and if we don't 
want to add another high energy scale to the model, 
the natural choice for compactification scale is M5 the 
fundamental quantum gravity scale which must be close 
to Electroweak scale ~ ITeV . It is easy to see that 
only if /i is of the order of Planck energy, the first KK- 
modes will not be produced in the present colliders Jlj 
or interaction of UHECRs in the atmosphere, otherwise 
KK-modes are light and apriori are easily produced in 
high energy collisions. 

In studying propagation of particles in the bulk we 
didn't consider the possibility and probability of pro- 
duction KK-modes in hadron collisions. Thus, ques- 
tion arises if the dominance of what is called small- a;^ 
events which are soft/semi-hard, permits the produc- 
tion of heavy KK-modes ^3 ■ For this purpose we need 
a model which can relate the low scale observables at 
small xb to the high energy scales where new physical 
phenomena such as escape to the bulk can happen. Our 
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argument about the need for such a model is that as we 
will explain in more details in the next section, at very 
high energies in majority of events, the signature of a 
new phenomenon is shrouded in low energy effects and 
we must have a suitable tool to extract the information 
from the mess of final states. 

In the following sections we first explain difficulties 
of analyzing the observations at very high energies and 
briefly review the QCD evolution models necessary for 
this purpose and in more details the new phenomeno- 
logical model called McLerran-Venugopalan Color Glass 
Condensate (MVCGC) model. Then, we extend this 
model to warped 4+1 space-times with brane bound- 
aries. And finally, we present the results for distribu- 
tion of gluons in the bulk. To be able to have a rough 
estimation about the propagation time in more complex 
brane models where the application of MVCGC model 
is more difficult, we also show some results from classical 
propagation of relativistic particles in the bulk. 

2. High Energy Physics 

Probing physics at very high energies is not an easy task. 
Even the clean signal of production of a heavy particle 
in the collision of two high energy leptons is dominated 
by elastic scatterings, production of light particles, con- 
fusions due to initial and final state soft radiative pro- 
cesses, etc. And we don't talk about technical difficulties 
of accelerating light leptons to very high energies ! It 
seems therefore that we are band to extract signals of 
new physics from billions of uninteresting interactions. 

The other possibilities - at least in the case of some 
models suggested as extension to the Standard Model 
- are the observation indirect effects and reconstruction 
of high energy processes by tracing back the low energy 
final particles to their high energy origin. The best ex- 
ample for the first case is the prediction of the mass 
of top quark from precise measurement of electroweak 
processes in LEP-1 with a center of mass energy much 
lower than the mass of top quark. The second possi- 
bility is more difficult and is less explored. Evidently in 
hadron colliders jets are reconstructed to obtain the rest 
mass of the partons which their hadronization make the 
jets and in this way it is possible to detect heavy parti- 
cles like top quark or Higgs. However, one can imagine 
other situations. For instance, the restoration of symme- 
tries broken at low energies can increase number of sym- 
metry charges - approximately massless gauge bosons 
with roughly the same coupling constant as low energy 
bosons. This is specially expected in presence of Super 
Symmetry when all the coupling become the same. To 
detect a signature of such a process, either one should 
select very rare events which permit direct observation 
of production of a heavy boson. Or one can try to find 
an evidence of its existence in events where such bosons 
are produced virtually or on shell and subsequently de- 
cay to secondary partons which their hadronization and 
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possible loss of a significant number of them out of the 
detector smears the signature of the heavy bosons. Dia- 
grams (O and are respectively examples of a dean 
and a messy: In the first diagram a heavy boson, 
labeled h , is produced in the collision of a pair of high 
energy collinear quark-antiquark and decays to a pair 
of quark-antiquark (or leptons). When it is close to its 
mass resonance, a significant enhancement in the num- 
ber of events with proper rest mass is observed. How- 
ever, the initial quarks (partons) must have enough en- 
ergy to reach the resonance. In the second diagram the 
production of the heavy boson - here assumed to interact 
with gluons after restoration of a symmetry e.g. the bro- 
ken symmetry responsible for SU{2) — SU{3) splitting 
- is smeared by feather decay to gluons which in turn 
produce other partons. The boson in this diagram can 
be virtual and therefore can influence the cross-section 
(evidently much less than the first diagram) even when 
the available energy is not enough to produce it on-shell. 

Even in the first diagram at very high energies the 
quark-antiquark from the decay of the heavy boson can 
in turn have QCD radiation and make low transverse 
energy jets which also interact with the remnants of the 
colliding hadrons and make the final signature of the 
heavy boson very weak. Therefor as the borders of un- 
known physics approach higher and higher energies, the 
probability of clean events becomes smaller and observ- 
ables are dominated by low energy physics. Note also 
that this depends on the nature of the extension of SM 
physics at very high energy scales and depending on the 
underlaying physics some effects are more difficult to ob- 
serve than others. It is therefore necessary to find a way 
to relate low energy observables to the hidden physics 
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behind them. 

At present, practically all the high energy collid- 
ers both man-made and cosmic colliders/interactions of 
high energy Cosmic Rays, are hadronic or semi-hadronic 
(lepton-hadron collision). It is well known that partons 
final state is always smeared by hadronization. Thus 
if we want to learn about the effect of exotic processes 
well before hadronization we must be able to relate the 
hadronic observables and final parton states to high en- 
ergy partons or other particles. This is not an easy 
task because of non-perturbative interactions and multi- 
particle state of both initial and final hadrons. The 
common practice is finding an evolution equation for 
partons. It should permit to trace back the observed 
jets in the final state as well as the rest of the rest of the 
hadron bag which does not participate in the hard in- 
teraction, to the point where their parent partons have 
~ 1/3 of hadron momentum. By approaching high en- 
ergy scales exotic phenomena such as appearance of an 
extra space-time dimension should appear as a devia- 
tion from Standard Model. As an example see |55] for 
evolution of gluon density in SUSY models. 

In hadron-hadron and hadron-lepton Deep Inelas- 
tic Scattering (DIS), the main observables which deter- 
mines the partial and total cross-sections are the den- 
sity distribution of partons and its evolution at differ- 
ent energy scale. It is well known that at medium and 
high colliding energies these interactions are dominated 
by a regime called small- a;^ where high transverse mo- 
mentum partons - mostly gluons - have only a very 
small fraction of the incident hadrons energy. The den- 
sity of gluons in this regime is very high such that, al- 
though the energy scale is much higher than Kqc d and 
cks < 1 , asln(l/a;B) > 1 and perturbative expansions 
fail. For the same reason parton density evolution equa- 
tions such as BFKL and DGLAP equations are not any- 
more valid approximations. The next subsection reviews 
very briefly the main aspects of these evolution models 
and their domain of validity. 

2.1. Partons Evolution Equations 

The non-perturbative and multi-body nature of QCD 
interactions in hadronic collisions makes the tracing of 
the final state particles to the parent partons very dif- 
ficult. From perturbative regime of the QCD, we know 
that structure functions - presenting the distribution of 
partons at each energy scale - have an approximate scal- 
ing properties. Their main variation is with respect to 
Bjorken parameter xb which can be interpreted as the 
fraction of initial hadron momentum carried by a par- 
ton. The QCD radiation corrections add a slow energy 
variation proportional to InQ^ (at lowest order) to the 
structure functions. Parameter is the absolute value 
of the square of invariant energy-momentum of the ex- 
change boson. More generally, it can be any energy 
describing the scale in which QCD interactions are hap- 
pening. When very high energy colliders were not yet 



available, the dependence on the transverse momentum 
fcj^ was usually ignored. In fact in DGLAP approxima- 
tion (see below) structure functions are integrated over 

faixB,Q^) = !k2^^Q2dkl<j)a{xB,kl,lnQ^) (6) 
a — gluon, quarks, antiquark. 

The reason for this mostly historical definition has been 
that at low energies and lowest order of perturbative 
QCD, k'j^ > A^Q^jj and without QCD radiation cor- 
rection, its maximum value is limited by . At high 
energies however, and fc^ become unrelated and the 
latter is controlled by the initial energy of the incident 
hadron. 

The core part of the structure functions at low ener- 
gies is non-perturbative and must be determined exper- 
imentally. Although due to non-perturbative effects the 
evolution with respect to In and fc^ can not be cal- 
culated from first principals of the QCD, good approx- 
imations can be found by using factorization of pertur- 
bative diagrams from non-perturbative regime and some 
statistical arguments. 

The most popular evolution equations (also by his- 
torical order) are DGLAP and BFKL. As we mentioned 
in the previous paragraph, DGLAP equation describes 
the evolution of structure functions with In : 

^j^^/a(a;i3,(5^) = -fa{xB,Q'^) 

[ dx xWa^b{x) + Y. I —Wb^a{x/y)My,Q^) 

I, Jo b y 

(7) 

The kernel Wa->b is the probability for parton a to 
change to parton 6 by a perturbative QCD interaction 
at a given xb and a given energy scale . The first and 
second terms in Q present respectively the probability 
of annihilation and creation of parton a in perturbative 
QCD interactions. 

At high energies and k\ are decoupled and the 
latter becomes a better indicator of QCD interaction 
scale. BFKL equation determines the evolution of unin- 
tegrated structure function (xs , fc^^ ) of gluons which 
are the dominant parton with respect to \n{l/xB)- Al- 
though BFKL equation is more suitable for determining 
the density of gluons at very small xb regime of high 
energy colliders, the calculation and interpretation of its 
kernel function is more subtle and the inclusion of higher 
QCD orders into its kernel is more difficult. 

Since HERA data at very small xb became avail- 
able, it has been realized that both DGLAP 26 and 
BFKL equations are inadequate for explaining the ob- 
served total cross-section and its evolution with kine- 
matic parameters. Moreover, both in HERA and in 
RHIC at very small x b , the cross-section which ac- 
cording to BFKL should exponentially increase, ap- 
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proaches a saturation and becomes much smaller than 
expected 

The main reason for the failing of DGLAP and 
BFKL is that none of them deal correctly with mo- 
mentum ordering In DGLAP the integration over 
transverse momentum is equivalent to considering only 
coUinear part and momentum ordering is completely 
washed out. For large xb the effect is negligible as the 
collinear component is much larger than transverse mo- 
mentum. But at small xb the ordering both in emitted 
and transferred partons becomes important. BFKL add 
the transverse momentum to the calculation of the ker- 
nels but assumes that all of them are of the same order 
(a regime called multi-Regge) j29| . 

Two strategies have been followed to find structure 
functions and evolution equations enough precise, spe- 
cially at small xb ,to explain observations. The first one 
called CCFM evolution equation is in the same spirit 
as DGLAP and BFKL models and uses the factorization 
of soft partons but takes into account the ordering of real 
and virtual momentums and leads to an evolution equa- 
tion which at large xb is correct to collinear approxima- 
tion (like DGLAP) and at small xb is not restricted to 
this approximation. Nonetheless, it has to make some 
simplifying assumptions about matrix elements in the 
kernel without which it is not possible to find a simple 
recurrent relation and thereby a simple evolution equa- 
tion. Parton distributions are unintegrated similar to 
(0 but also depend explicitly on . The kernel (split- 
ting function) includes form factors for IR regularization 
of soft on-shell and virtual giuons, latter is particularly 
important at low x b ■ The differential representation of 
this evolution equation is the following [221 j28| : 

q'^d A{xB,k±,q; ko±,Q, 



^ dz P{z, ^,k±) XB q 



z A,(f,Q) 



A(^,k'^,l;ko^,Q,^^) (8) 



where k'j_ ^ kj_ + ^q. The splitting function of CCFM 
evolution equation P is complex and depends on the 
kinematic of emitted partons: 



P{z,q,k±) = JJPi(zj,gi,A:i_L) 



(9) 



^s{kf^ 



A„,(z„g„A:,^) + "^^^^ ^'^''^'^ (10) 
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^ns{zi, qi, ki±) = 

dz' d(f 

2 ^ 



1 - z 



exp -as 



e(fc,i - q)e{q - z'qi) 



(12) 



The second method is quite different in spirit and is 
based on a phenomenological modeling of many-particle 
QCD processes. The initial idea raised by L. McLerran 
and R. Venugopalan [301 in the context of determina- 
tion of partons density in a large nuclei and is based on 
a classical treatment of color fields. Later however it 
was observed that it can be also applied to other situa- 
tions such as small x b regime in collision of high energy 
partons where the density of giuons is very large and al- 
though in perturbative regime of QCD, large probability 
of interaction make the non-perturbative effects impor- 
tant. For this reason this model is now called Color 
Glass Condensate ( CCC). The reason for words Color 
and Condensate in this naming is clear. The word Class 
is used because giuons at small xb are produced from 
giuons with high xb i.e. ones with a larger share of 
the energy and momentum of the initial hadron. In infi- 
nite momentum frame they have a time dilatation which 
will be transferred to lower energy scale partons. They 
evolve slowly compared to natural time scale, similar 
to a glass behavior [3Tj . This observation is very im- 
portant for extending this model to higher dimension 
space-times which should manifest only at very high en- 
ergy scales. 

It is important to mention that the condensation of 
giuons in this model is only the asymptotic state of gluon 
matter at low energy scales and high densities, and the 
model can be well used at higher energy scales where 
the density of giuons is not enough high to make a con- 
densate. This can be seen from the possibility of finding 
BFKL approximation at the lowest order in this model. 

The idea of this model comes from the observation 
that a hadron with large momentum in the rest frame of 
the observer is contracted in the direction of the boost 
and its content is concentrated on a sheet-like surface. 
li E ^ Aqcd , partons are freed from confinement and 
can move in the direction of the boost. Due to interac- 
tion with the color sheet and themselves, their lifetime 
is however very short with respect to the time varia- 
tion of the fast partons on the sheet. When two such 
hadrons with opposite momentums collide with each 
other, just before the collision most of the swarm of 
slow partons (mostly giuons) is concentrated between 
two sheets made by fast color charges of each hadron 
(see Fignj. In a classical field theoretical view, one can 
assume the sheets as an external charge configuration 
playing the role of a source for the QCD gauge field 
(gluon field) propagating between two hadron sheets. 
This is a QCD analog of a capacitor in electromag- 
netism. One can therefore solve the field equation and 
determine the field strength (gluon density) for this con- 
figuration. The main difference between this model and 
its electromagnetic analog is that in the latter case the 
charge density of the sheets is well defined. In QCD the 
charge depends on the energy scale, in other word de- 
pends on how close or far the observer is looking at one 
of the sheets (hadrons) because relating a charge to the 
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Figure 1: Collision of two hadrons. The color /grey level 
of the dots is an illustration of color charge of the gluons. 
The thickness of sheets presents the scale A+ for which 
the color charge is defined. 



sheet or to the swarm is not unique; closer an observer 
looks at the hadron, more of gluons s/he is seeing be- 
long to the swarm in the space between sheets and less 
to the charge sheets 

This phenomenological illustration is very helpful 
when we want to extend hadron-hadron collision to the 
very high energy scales where exotic processes such as 
extra-dimensions can be accessible to partons [32] ■ For 
instance, if hadrons are boosted such that the physi- 
cal width of the hadron sheets becomes comparable to 
the effective size of the bulk, at very high energy scales 
when slow partons leave the charge sheet, they can enter 
to the bulk - assuming that at these scales there is no 
difference (symmetry breaking) which makes the visible 
brane the prefered direction. According to Glass con- 
cept explained before, the partons at high energies are 
the parents of low energy, small x b partons and thus by 
escaping to the bulk naturally they take with them the 
low energy partons. This leads to a significant reduction 
of number of partons on the brane which must be ob- 
servable. Even when gravitational force of the warping 
can bring back the particles to the brane, a time deco- 
herence - equivalent to a small mass excess should be 
observable. Notice that this picture is Lorentz invari- 
ant. To an observer in the frame of one of the hadrons, 
partons of the other hadron come from very short dis- 
tances and can participate in exotic processes including 
emission to extra-dimension. 

Although it is possible to add the effect of extra- 
dimensions with non-flat metrics to parton evolution 
equation such as BFKL and CCFM, we have prefered 
to use the CGC model for this purpose. The reason 
is the consistent structure of this model which, as we 
will show in the next section, is manifestly Lorentz in- 
variant. Moreover, at least formally, there is no need 



for many simpliflcation assumptions in the intermediate 
steps which is the case for all the evolution equation 
formalisms. 

3. Color Glass Condensate (CGC) in 
4+1 Space-time with Brane 
Boundaries 

In this section we use light-cone coordinates defined as: 

x+ =x^ + , X- ^x° -x^ (13) 

In McLcrran-Venugopalan approximation, QCD inter- 
actions are modeled by an effective SU(2i) gauge field 
of gluons. For simplicity we neglect the color index 
except when its explicit indication is necessary. In Light 
Cone (LC) gauge = 0. The time variation scale of 
color charge density on the sheet produced by the in- 
coming hadron is much longer than the gluon swarm 
time scale. Therefore, it can be considered as a static 
charge. The classical dynamic equation of the effective 
gluon field is: 

[Da,F^^]{x) = S^+W{x+,x)p+{x)W\x+,x)(U) 
W{x+,x) = Texplig dri+ — {r]+ , x) I 

(15) 

Da = d,A-i9A'XT- (16) 

The symbol ; means covariant derivation with respect 
to RS metric. Color charge p is -I- component of color 
current: 

J^{x)^S^+p{x-,x^,z) (17) 

The Wilson-line term in the right hand side of (|14|l guar- 
antees the gauge invariance of this equation. The time 
ordering operator T operates in (|15|l on A^ = A^T"^ 
and orders fields from right to left in increasing sequence 
of x+ . 

The subscript A+ of the color charge density means 
that it is defined at scale A"*" — xbP'^ where P"*" is the 
energy of the initial hadron. The scale A"*" indicates the 
maximum LC momentum of the parton swarm. In clas- 
sical MV model it is not possible to relate models at dif- 
ferent scales. But when quantum corrections are added, 
one can determine the evolution of parton density with 
scale A"*" and thereby relate the gluon distribution at 

— > to xb — > 1 and vice versa. Physical quantities 
such as 2-point correlation functions at classical (tree) 
level are obtained from solutions of (|14|l integrated over 
all possible p{x) with probability distribution Wa+ [p] ■ 

We assume a universal coupling in the bulk and on 
the branes. The aim of having additional coupling on 
the branes in |15| is the request for confinement of vec- 
tor fields on the brane when fermions are confined. It is 
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however easy to see that in this model fermions can't be 
confined to the visible brane because this violates the 
gauge invariance and makes the model inconsistent. 

In LC gauge equation H14|) has a solution with — 
in addition to the LC gauge condition yl+ ~ |33|. 
The solution for other components is: 



Ai — -U{x)d-iUHx) i^x^,z 
9 



(18) 



with Z// an element of Q CD SU{3) group. Unfortunately 
in this gauge there is no analytical way to relate U to 
the charge p and to find out its explicit form. But in 
covariant gauge one can find an explicit solution. If 
we perform a gauge rotation to covariant gauge, A^ 
becomes: 



A^ = U^A^U + -U'^d^U 



(19) 



In LC gauge A'^ = and has Kahler potential form: 



A+ = -U^d+U 
9 

and equation reduces to 

-a,d^A+ = p{x) 

p{x) = U\x)p{x)U{x) 



(20) 



(21) 
(22) 



p is color charge density in covariant gauge. Using (|^ . 
equation (j2()|l can be inverted and: 

L({x~ , x^ , z) = P exp< ig / dry^— ^+(77^, x^, z) 
L Jxg ^ 

(23) 

where P orders A~^ = A^T°- from right to left in in- 
creasing or decreasing order of x~ argument respec- 
tively for x" > Xq or x^ < ■ We need both LC and 
covariant gauge formulations because in the latter gauge 
the solution of 1)14(1 is simpler, but the gluon distribution 
function has a simpler description in LC gauge. 

Solution of gives the propagator (Green func- 
tion) in the warped 4-1-1 space-time. After defining the 
Fourier transform of A~^ : 



i+(fc^,z) 



(24) 



one can expands equation ((21|l and obtain the equation 
for the propagator: 

z^d^,A{z, z')+zd,A{z, z')+(fc'z'-l)A(z, z') = RH{z~z') 

(25) 

with k'^ = — fc^ . General solution of (|25|) is: 

A(z, z') = C{z')Ji{kz) + D{z')Ni{kz) (26) 



where Ji and iVi are respectively Bessel function of first 
and second type. Note that the momentum of partons 
in the direction of the initial hadron i.e. A"*" = xbP^ 
appears only in the scale of the model or equivalently 
in the color charge p. Also the dependence on is 
not dynamical which means that it should be fixed when 
boundary conditions are imposed. Integration constants 
C{z') and D(z') are determined by applying boundary 
conditions on the branes at z = i? and z — R' . This 
results to following mass spectrum for KK-modes for 
large n: 



R' 



L ^ 37r 



(27) 



Therefore the coefficient Xn in CJ for lightest modes 
is ~ 2.36. The mass of all the gluon KK-modes in this 
model is real. The spectrum begins from a massless zero- 
mode and there is a gap between zero-mode and higher 
modes proportional to p! = R'~^ which for all macro- 
scopic bulk is very small. In fact for large pL > 30 , even 
for compactification scale p ~ Mpi , the mass of lightest 
KK-modes is much smaller than CM energy of UHECR 
interaction in the atmosphere and when the interaction 
scale A+ > m„ , KK-modes can be produced. 

There is also a zero-mode with fcj^ = with A(z, z') 
satisfying: 



d2 

zd,{zd,A{z, z')) - 1 = ^5{z ~ z') 



A{z,z') = Co{z')z + Do{z')z-^ 



(28) 



(29) 



After applying the boundary conditions to 129() . one 
finds that the integration constants Co(z') and Do{z') 
are in general non-trivial and therefore the zero-mode 
propagates to the bulk but its wave function exponen- 
tially decreases inside the bulk. This is a well known 
fact nil and contrasts with scalar field case in which 
the Green function of the zero- mode is cc S{z — z') and 
consequently induces a discontinuity in the spectrum 
and separates zero-mode from higher KK-modes. 

Propagation of gluon swarm to the bulk in MVCGC 
model has important implication for attempts to local- 
ize massless gauge bosons by adding an induced kinetic 
term on the branes. It has been shown [221 El that 
this term increases the coupling of the zero-mode to it- 
self and if fermions are confined to the brane, the gauge 
vector field becomes quasi localized. Here it is easy 
to see that due to z dependence of and therefore 
U{x) , according to H22|l . p will depend on z even when 
p is confined to the brane. We will see later that be- 
cause of the renormalization group equation which re- 
lates Wa+ [p] the distribution of color charge density at 
each A+ scale to other scales, p or p gradually receives 
a z dependence, i.e. the color charge escape to the bulk 
even if initially - at large A+ - it is concentrated on the 
brane. A more physical image of this process is obtained 
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if one considers the transverse momentum ordering. At 
smaller 's partons have larger transverse momentum 
because they have lost their energy by QCD radiation 
which at the same time increases their transverse mo- 
mentum. When an extra space dimension is available, 
with each radiation the transverse momentum of the 
outing parton in the extra-dimension direction increases. 

Finally gluon distribution can be calculated from the 
matrix elements {A\^A^i^) after quantizing the classical 
solution A\ ■ The definition of the gluon distribution 
function (structure function) in LC gauge |34| can be 
extended to 4-t-l space-times and results to the following 
expression: 



xbG{xb, Q', z) = 2Ji{1)P+xbQ'' j dP 



„-L2 



dy'e 



'e^y' (A\x+,x',y')A{x+,x-x',y-y')^ 



(30) 



In this definition of gluon distribution function is the 
3-dim transverse momentum scale and xbG{xb,Q^, z) 
is the density of soft /semi-soft gluons up to \k±\ < Q . 
Due to the bounded non-flat geometry of the bulk we 
can not extend \k^_\ to the whole transverse directions. 
The brackets () in H3()|l present both quantum mechani- 
cal expectation function and averaging over all possible 
configuration of color charge density p{x) . The proba- 
bility functional Wj^+ [p{x)] for each configuration at a 
given scale A"*" can be obtain from path integral quan- 
tization of MVCGC model. It has been proved that 
Wj^+ [p{x)] evolution satisfies a renormalization group 
equation P5): 

SWr[p] (1 



5t 



2 5p'^{x^,z)5p^^{x'^,z 
5 



6p'^(x^, z) 



(31) 



where r = ln(P+/A+). Indexes a and h are color in- 
dexes. Matrix elements ct" and x"** are defined as: 



X 



a'^ = {Sp^ 
" = {Sp^Sp') 



(32) 
(33) 

where is the fluctuation of p distribution from its 
classical value. Lowest order QCD diagrams which con- 
tribute to charge variation are discussed in |35j . 

In Gaussian approximation (see (|38|l below) the 
standard deviation (t^+ ~ ^afcX"'' at lowest order 
(BFKL approximation) |35|: 



-g'D+Tr T^Gl,,ix,y) 



(34) 



x=y 



'-'0 



Xab 

\x) 



A^g^T+^{x){x\Go^M^+\y) (35) 



i — z 



(36) 



where Gq is the propagator of free gluons. In equation 
(|35|) is the gluon classical field strength of solution 



of H14|l . This approximate evaluation of x can be used 
in (|37|l below to determine the gluon distribution up 
to BFKL approximation. It can be shown 32 that in 
this approximation the matrix elements {Al^AD have 
the following expression: 



'^aC'^' z)Abi{x , Z ^ 



SabXix, X', z, z')d!d'.aix^,x'^,z, z') + ... (37) 



where: 



R 



(Ta+ [x ,X ,Z) 



dx" a\+{x" ,x'^,z) 
(38) 



and: 



{d,kd:'')-''S'\x^ ~x'^,z,z')=j{x^,x'^,z,z') 
1 



(27r)2 
^' dz" 

^A{z",z,k^)A{z",z',k^)+j{0) (39) 

R ^ 

The standard deviation of the Gaussian distribution 
f7A+ can be replaced by BFKL approximation (|32|l . 

The result of numerical calculation of the gluon 
distribution for the fine-tuned RS model and when 
the bulk scale p is of the same order as fundamen- 
tal gravity scale M5 are respectively shown in Fig|5] 
and Fig.fig:rshighmu. They shows that in both cases 
the probability of gluon emission into the bulk is much 
larger than their emission into the brane. Approxi- 
mations we made to obtain these distributions have a 
validity domain roughly the same as BFKL approxima- 
tion. Nonetheless, they can be evolved to lower scales 
by using (|31|l . Even without any additional calculation 
and just by remembering that gluons at lower scales are 
produced from gluons at higher scales, one can conclude 
that most of the softer gluons should also be emitted in 
the bulk. Thus, if il/5 < lO^^eV^, one had to observe 
a much lower flux for the high energy wing of UHE- 
CRs. In fact not only it does not seem to be the case, 
the observed flux is much higher than classical expecta- 
tions and either a top-down model of a decaying heavy 
dark matter or extremely exotic astronomical sources 
are needed to explain their flux. 

4. Classical Propagation of Relativistic 
Particles in the Bulk 

Now that we have shown that high energy particles in a 
hadron-hadron collision penetrate to the bulk, we would 
like to use the classical propagation of these particles in 
the bulk to constrain other brane model. The reason 
is that although the MVCGC formalism or other evo- 
lution models such as CCFM can be extended to any 
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D. 10 20 30 5 10 15 20 
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Figure 2: Gluon distribution in the bulk normalized to 
the amplitude of the distribution on the visible brane at 
z = R' for Randall-Sundrum models for 2 x lO^ey ^ 
\k±] ^ 1.26 X IQ^^eV . Curves are rainbow color/grey 
coded. From bottom to top M5 = lO^^eF, lO^eF and 
lO^^ey , with fi obtained from fine-tuned RS model (O 
i.e. /i = lO^^^eV, \Q^^'^eV and \Q~^^eV respectively. 
For left plots log(i?'/-R) = log(Mj,//Af5) , i.e. the same 
as (@J). In right panel a smaller R' / R is used to see the 
effect of a smaller bulk with same fundamental gravi- 
tation and compactification scale. Our tests show that 
apparent vibration of the distribution is mainly due to 
low resolution of our numerical calculation. 



brane model, for more complex metrics the calculation 
becomes very difficult. For a simple order of magnitude 
estimation of parameter space, a time of flight estima- 
tion can be adequate. The general form of the metric of 
a 4+1 dimensional space-time with curved bulk is: 



n^{t, y)df-a^{t, y)5,jdx'dx^ -b^{t, y)dy^. (40) 



In a gauge where b{t, y) = 1, after solving Einstein 
equations and applying boundary conditions on the 
branes, one obtains the following solutions for a{t, y) 
and n(t, y) : 

a^{t, y) = A{t) cosh{fj.y) + B[t) sinh(^y) + C{t) (41) 









Figure 3: The same as Fig|5|but with = Afs. Other 
details are the same as Fig|5] In the top right plot 
the difference between curves for different k-^ is much 
smaller than the resolution of this figure and they are 
overlapped. 



a'{t,y) = n''{t,y)allt) = 



A{t) cosh(/xy) + B{t) sinh(^?/) + C{t) 



4a2(t,y) 



A{t) 
Bit) 

C{t) 



n{t,y) 



aQ^t)-Cit) 

-Pbo«o^(i) 
2dlit) 



2^2 

a{t,y) 

doit) 



Pb\ 



(42) 
(43) 
(44) 

(45) 

(46) 
(47) 



For any density p, p' = p/Ans, Ars = •ip./k^- The 
densities and ps are respectively effective total en- 
ergy density of the brane at y = and the bulk. We 
consider only AdS bulk models with pB < 0. The con- 
stant = 87r/M| is the gravitational coupling in the 
5-dim. space-time. The model dependent details such 
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as how pI,^ and ps are related to the field contents in 
the bulk and on the brane and how they evolve are irrel- 
evant for us as long as we assume a quasi-static model. 
The solution (|41|) is valid both for one brane and multi- 
brane models. The only difference between them is in 
the application of Israel junction conditions [201 |37| . 



After changing the variable y to z 
(133: 



e^"^ and using 



V 



2z 



V{t,z) 



(48) 



2C{t) 



z + 



„2 / 



dz 



(49) 



(50) 



9 is an integration constant. If an ejected particle to 
the bulk comes back to the brane, its velocity in the 
bulk must approach to zero at its turning point before 
the particle arrives to the bulk horizon (if it is present). 
The roots of (|50l) correspond to these turning points and 
determine the propagation time in the bulk. The typical 
propagation time we are interested in is much shorter 
than the age of the Universe and therefore B, C and 
A^B^C during propagation are roughly constant. The 
right hand side of ()50|l depends only on z and is easily 
integrable: 



2(istop~io) 



2dz 



20 



(51) 



In H51|l . to is the initial time of propagation in the extra- 
dimension and tstop is the time when the particle's ve- 
locity changes its direction, i.e. when dz/dt ~ 0. The 
integral in H51|) is related to the elliptical integrals of the 
first type J-{LiJ, v) where oj and v are analytical func- 
tions of the denominator roots in H51() and zq . Note 
that Zstop corresponds to the closest root to zg . 

Using this simple model, we determine the propaga- 
tion time of relativistic particles in the bulk for some of 
the most popular brane models. This time delay must 
be smaller than the time resolution of air-shower detec- 
tors as no time decoherence - delay between arrival of 
most energetic particles - has been observed. 

FigElshows Atpropag 8ls a function of /iL and Afs for 
massive relativistic particles. With present air shower 
detectors time resolution of order 10~^ sec, only when 
M5 > lO^^eT^, the model is compatible with the ob- 
served time coherence of the UHE showers. For fine- 
tuned RS model n « G/k^ ^ i.e. fx = GM^ ~ IQ-'^eV 
for M5 ~ lO^^el^ 13 ■ Due to smallness of /x and con- 
sequently lightness of KK-modes for SM particles even 
this model with large M5 has already been ruled out 
unless a conserved quantum number prevents the pro- 
duction of KK-modes [HH|- For massless particles: 



z{t) - zo = p.{t - tof 



(52) 




100 



Figure 4: Propagation time for relativistic particles 
with Lorentz factor u"(^o)/^ — 10^ (full fine) and 
u°^{to)/N = 1.2 (dash line) in RS model. From 
top to buttom curves correspond to M5 = lO^^eV , 
Ah = IQi^ey and M5 = W^eV (or ^ ~ IQ-^el/, 
/i ~ 10~^^ey and /j, ^ 10~^eV for fine-tuned model) 
respectively. The horizontal line shows the time coher- 
ence precision of present Air Shower detectors. 



In (|52|l . z{t) is monotonically increasing and there is no 
stopping point. With our approximations there is no 
horizon in the bulk because a(t, y) is roughly constant. 
Therefore, equation ((^21) means that massless particles 
simply continue their path to the hidden brane and their 
faith depends on what happen to them there. At very 
high CM energy of UHECR interactions if charged par- 
ticles can escape to the bulk photons are also dragged 
to the bulk and never come back. This is very similar 
to the conclusion of more precise and completely inde- 
pendent formulation of the previous section and shows 
that the logic behind them are consistent. 

Next we consider the case of a general 2-brane model. 
Numerical solution of constrained 2-brane models in 
|SZ| shows that for fiL > 5 the tension on both 
branes can be positive and very close to Afis ■ We use 
constraints on the Cosmological Constant and hierar- 
chy to find Pq and p'^ = p'^ , the branes tension. We 
redefine them as p'q = 1 + App and p^ = 1 + Ap^ . To 
solve hierarchy problem the following constraint must 
be fulfilled: 



M| 



n\ PAo(1 ^ cosh(pL)) + sinh(pL) 
Uq p'k (1 — cosh(pL)) -I- sinh(pL) 



< 1 



(53) 



This leads to: 



Ap^ = 



iV^ 1 _ g-''-^ + Api^(l - cosh(pL) 



1 — cosh(pL) 



1 - g-*'-^ 
1 — cosh(pL) 



(54) 
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For a very small N'^ and Ap^ < 1 , the first term in 
(El is 0{N'^) and: 



1 



1 - e-^'^ 



1 — cosh{fiL) 1 — cosh(/iL) 



2e'f^ (55) 



Using a\/a\ — where H is the Hubble Constant 
on the visible brane |37j: 



1 



27V2 sinh(/iL) 
(1 — cosh(^L)) 



- 1 ± 



(1 — cosh{^L)) 



2H^ 



+ e 



— fiL 



sinh(^L)(^ + 2e-^^) 



In H56|l the solution with plus sign gives Ap^ 



(56) 



-2 



which deviates from our first assumption |Ap^| < 1 and 
leads to a negative tension on the visible brane like static 
RS model. The solution with negative sign is Ap'^ w 
2e~'''^ and both branes have positive tension close to 

Using these values for the tensions in (|41|) to (|47|l . 
we can find the time delay for this class of brane models: 



Atr, 



M8|Ap^ + C'|) 



2 arctan / 



/ g(Z- - Zq) 

p{zo - z+) 



(57) 



(58) 



1/. 2 

= - 2 H 

2\ pq 

C'{C' - 02) + 4(2 + Apf, + C'){Ap'a + C 



2 _ 



c 



c 



2 + Ap'a + C 
Ap'o + C 



-C = 



Ap'o + -{N' 



(59) 
(60) 

(61) 

(62) 
(63) 



In H58|l . Zq = e^^ and, z± are the roots of velocity 
amplitude. The result of the numerical calculation of 
the time delay for these models is shown in FigEl We 
conclude that for these models the time delay can not 
constraint the value of Afs , but put constraint on p,L , 
roughly speaking the warp factor. This means that by 
adjusting the warping, it is always possible to make the 
time delay too short to be observable at present air- 
shower detectors. Nevertheless, in these models fi is 





1 00 1 50 200 



100 150 200 

mu L 



Figure 5: Left: Propagation time for relativistic parti- 
cles in 2-brane model. Description of the curves is the 
same as Fig0] Right: Parameter p{eV) as a function 
of piL . It is roughly independent of M5 . 



not arbitrary and depends on pL . On the other hand, 
in many radion models for stabilizing the bulk, p is re- 
lated to the mass of radion and its value must be in the 
range requested by the models for stabilization. There- 
fore comparing the results presented here permits to in- 
vestigate which type of stabilization and particle physics 
behind them can satisfy all the constraints. Test of more 
models can be found in 

5. conclusion 

In this proceedings we tried to show that physics in fu- 
ture colliders and air-shower detectors is crowded by low 
energy effects and great efforts must be dedicated to ex- 
tract information about new physics hidden deep inside 
collision events. 

We investigated the effect of the escape of partons 
to a warped bulk. We showed that it can significantly 
reduces the total cross-section and multiplicity - number 
of final state particles emitted into the brane. To arrive 
to this conclusion we needed to use an evolution equa- 
tion to related the low energy observables to high energy 
parton (gluon) distribution in hadron-hadron collisions. 
We used the recently developed phenomenological model 
Color Glass Condensate and extended it to 4-1-1 warped 
space-times with brane boundaries and static Randall- 
Sundrum metrics. 

Inspired by the observation that relativistic particles 
can be emitted to the bulk, we have also used a simple 
classical model of propagation of relativistic particles in 
bulk for more general type of brane models and obtained 
constraints both for the fundamental scale of gravity and 
for parameters of the brane models. 
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